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This paper studies the finite temperature Casimir force acting on a rectangular piston associated 
with a massless fractional Klein-Gordon field at finite temperature. Dirichlet boundary conditions 
are imposed on the walls of a d-dimensional rectangular cavity, and a fractional Neumann condition 
is imposed on the piston that moves freely inside the cavity. The fractional Neumann condition gives 
an interpolation between the Dirichlet and Neumann conditions, where the Casimir force is known 
to be always attractive and always repulsive respectively. For the fractional Neumann boundary 
condition, the attractive or repulsive nature of the Casimir force is governed by the fractional order 
which takes values from zero (Dirichlet) to one (Neumann). When the fractional order is larger than 
1/2, the Casimir force is always repulsive. For some fractional orders that are less than but close to 
1/2, it is shown that the Casimir force can be either attractive or repulsive depending on the aspect 
ratio of the cavity and the temperature. 
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I. INTRODUCTION 



During the past decade, fractional calculus has attracted considerable attention from physicists and engineers. 
^ ■ This is particularly true for researchers working in the fields of condensed matter physics, where fractional differential 
[ equations have been used to model various anomalous transport phenomena Jj, S S) IE] ■ Thus one expects fractional 
calculus, in particular fractional differential equations, to play an important role in quantum theories of mesoscopic 
^ ' systems and soft condensed matter which exhibit fractal character. Fractional dynamics provides a natural framework 
^ ■ for describing the evolution of physical systems in fractal and multifractal media. By extending the above argument to 
■ quantum theories in fractal space-time, one may then have to deal with quantum mechanics and quantum field theory 
which satisfy fractional generalizations of Schrodinger, Klein-Gordon and Dirac equations. However, applications of 
fractional calculus to quantum theory are still relatively new. In quantum mechanics, fractional Schrodinger equation 
[ in its various forms has been studied recently 0, d, S [Tot . Fractional supersymmetric quantum mechanics has 
\^ • also been considered [llj. Fractional Klein-Gordon equation [ll|, [l^, [H, [13, flslTlGt and fractional Dirac equation 
' [13, m, [il] have been considered more than a decade ago for various reasons. For examples, field theory with 
0\ La grangian containing nonlocal kinetic terms involving fractional power of D'Alembertian operator arises in the (2 
-|- l)-dimensional bosonization 



in the underlying local theory 



20 J |21| and also in effective field theory with some degrees of freedom integrated out 
23 . l23j . Some studies have indicated that the QED radiative correction leads to a 
k>( modification of the propagator of a charged Dirac particle which acquires a fractional exponent connected with the 
rS fine structure constant that is a fractal propagator [24. . ,25 , 26] . In quantum field theory of gravity with asymptotic 
■ safety j27|, spacetime geometry cannot be understood in terms of a single metric. There is a need to introduce a 
different effective metric at each momentum scale. In addition, spacetime structure in asymptotically safe quantum 
Einstein gravity at sub-Planckian distances has been shown to be fractal [1^, |2^ . Numerical simulations have studied 
the propagation of a scalar particle in a dynamically triangulated spacetime, with a discretized version of the Einstein 
dynamics, and found that the spectral dimension of the microscopic spacetime is two and it tends to four for long 
time scales [sO, [111, [S^l . The use of differential and integral operators of fractional constant and variable orders may 
seem appropriate and natural for such quantum field theories in fractal and multifractal spacetime. 

Recently, work on Casimir energy associated with fractional Klein-Gordon field has been carried out, in particular 
the Casimir effect associated with fractional massive and massless fields has been studied [s^ [sj] ■ Self-interacting 
scalar massive and massless fractional Klein-Gordon fields have been considered in the context of topological symmetry 
breaking on toroidal spacetime (35.] . In this paper, we want to study Casimir effect corresponding to fractional Klein- 
Gordon massless scalar field in the piston setting. Casimir effect associated with piston geometry has attracted 
considerable attention since it was first introduced by Cavalcanti [s^. The geometric setup of Casimir piston can be 
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FIG. 1: A three dimensional piston. 



used to avoid divergence problems that usually plague the calculations of Casimir effect. By taking suitable limits, 
the piston approach can be used to derive the Casimir force acting on two parallel plates embedded orthogonally 
inside an infinitely long chamber [STj , and the Casimir force acting on two infinite parallel plates (ssl . [s^ . Therefore 
it would be interesting to study the behavior of the Casimir force acting on a d-dimensional rectangular piston due 
to fractional Klein-Gordon massless field. In line with our consideration with the fractional field, we also impose on 
the piston fractional Neumann boundary condition which allows the interpolation between the ordinary Neumann 
boundary condition and the Dirichlet boundary condition. 



II. CASIMIR PISTON FOR FRACTIONAL KLEIN GORDON FIELD WITH FRACTIONAL 

NEUMANN CONDITIONS 

Let ^(x, t), X e M'^, i e M be a massless fractional Klein-Gordon field with Lagrangian 

£ = i0(x,t)(-A)"(/.(x,i), 

where A = ^ + X^^Li ^pr- For the interpretation of the fractional operator (—A)", we refer to our previous work 
[33 |. The purpose of this paper is to investigate the finite temperature Casimir force acting on a piston moving 
freely inside a rectangular cavity [0, Li] x ... x [0, L^] (see Figure [T]) due to the fractional field 0(x, i). We assume 
that the piston has negligible thickness and is orthogonal to the xi direction, with position given by xi — a, where 
< a < ii. On the walls of the rectangular cavity, the field (/)(x, t) is assumed to satisfy Dirichlet boundary conditions, 
i.e. (?!)(x, t)\x =o or L- = 0, 1 = 1, . . . , d. On the piston, the field 0(x, t) is assumed to satisfy the fractional Neumann 
condition 

= 0. (1) 

Xi—a 

The cases where 77 — (Dirichlet boundary condition) and t] — I (Neumann boundary condition) have been discussed 
in our previous work [40|. Therefore in this paper, we only consider the case where 77 G (0, 1). 

Assume that the system is maintained in thermal equilibrium at temperature T. The Casimir free energy of the 
system is the sum of the Casimir free energy Eq^^ inside Region I [0, a] x [0,^2] x . . . x [0, L^], the Casimir free 
energy -Eq^s inside Region II [a,Li] x [0, L2] x ... x [0, Ld], and the Casimir free energy -Ega* outside the cavity 
[0,Li] X . . . X [0,Ld], i.e., 

ppiston _ 771I , 771II , ™ut 
-'^Cas - ^Cas + ^Cas + -C^Cas- 

The Casimir free energy of the exterior region does not give rise to Casimir force on the piston. Using Matsubara 
formalism, the Casimir energy inside the cavity [0, a] x [0, L2] x . . . x [0, Ld] is given by 

T 

£^Cas(a) = ^logdet 

where /i is a normalization constant with dimension length^^, det([— A//!'^]") is the determinant of the fractional Klein- 
Gordon operator [— A//i^]" acting on functions 0(x, t) on [0, a] x [0, L2] x . . . x [0, Ld] x [0, l/T], which satisfy Dirichlet 



dx^' 



[Xl,X2, 



,Xd,t) 
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boundary conditions on the boundaries Xi = 0, 1 < i < d, Xi = Li, 2 < i < d, satisfy fractional Neumann condition 
(III) on the bounclary xi = a and satisfy periodic boundary condition in the t direction, i.e., (/) (x, i + =0 (x, t) . As 
is discussed in jsj] , a complete set of eigenfunctions is given by 



(l)k.i{x,t) = sm j [[sm ^— 



'•^J 1 2iTilTt 



with eigenvalues Luj, ^//x 



, where 



^1,1 = 

\ / i=2 ^ / 

Here k = (fci, /c2, ■ • • i fed); fci € Z, ^2, . . . , fed G and Z e Z. The Casimir energy £'cas('^) t'^^^ given by 



We regularize this sum by exponentially cut-off method, i.e., 

i^k.(«;A) = f EE(i°g^)^^'"'"''' (2) 

where A is a cut-off parameter with the dimension of length. For the Casimir energy in Region II, one can show that 
it can be obtained from the Casimir energy in Region I ([2]) by replacing a with Li — a. Therefore, using the result 
(jA3P in the Appendix, we find that up to the term constant in A, 

Eh,,{a; A) + E'^,,{a; A) = -aT^ ^c,+i_.(ii) (log[AA.]2 - 2^(z)) A"^ - ^ (c'(0; a) + C'(0; Li - a)) , (3) 

»=2 Mij ^ ^ 



where 



2<ji<j2<...<id-i-.<ti 



and ^(s; a) is the zeta function 



C(s;a) = E E '^fc.? 



2s 



The first summation in ^ is the sum of the A ^ 0+ divergent terms. We notice that these divergent terms are 
independent of the position of the piston a and the fractional parameter rj. The last term in ([3]) is what one would 
obtain for the Casimir energy if one uses the zeta regularization method to compute the Casimir energy. 

Upon differentiating ([3]) with respect to a, we find that all the terms that diverge when A — > 0^ do not contribute 
to the Casimir force acting on the piston, since these terms are independent of a. By taking the limit A — s- 0"*", we 
find that the Casimir force acting on the piston is 

Fcas(a; L,)=- ^hm |^{i?Cas(«; A) + SlUa; A)} = (C'(0; a) + C'(0; L, - a)) . (4) 

Using the result (|A4|) in the Appendix, we find that the Casimir force can be written as the difference 

Fcas(a; Li) = Fgi^ia) - Fgi,{L, - a), 

where 



'Et2(ir) +(27rff)2 

FSLsi^) = -2"™e E E / / ; X (5) 
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is the Casimir force acting on the piston when the right end of the cavity is infinite distance away, i.e., Li — > oo. 
It can also be interpreted as the Casimir force acting between two parallel plates embedded inside an infinitely long 
rectangular cylinder, where Dirichlet boundary condition is imposed on one of the plates, and fractional Neumann 
boundary condition is imposed on the other. It should be remarked that by putting 77 = and 77 = 1 in ([5]), we obtain 
respectively twice the Casimir force when Dirichlet and Neumann boundary conditions are imposed on the piston. 
The factor of two arises since for Dirichlet and Neumann boundary conditions, the ki and —ki modes are dependent, 
but for general fractional Neumann conditions with 77 e (0, 1), the ki and ~ki modes are independent. 
Since 



Re 



X cos(7r?7) — 1 



(x — cos(7r77))^ + sin (irri) ' 



(6) 



we can conclude from (O that the Casimir force is always repulsive when 77 > 1/2. For rj e (0, 1/2), the Casimir force 
can be attractive or repulsive depending on the values of a, L2, . . . , Ld and T. However, when a is large enough, or 
more precisely, when 



a > 



log cos(7r?7) 



the Casimir force is always attractive. Taking the derivative of f{x) ([6]) with respect to x, we have 

— a;^ cos(7r?7) + 2x — cos{TT'ri) 



((x — cos(7r7/))2 



sin^(7r77))^ 



From this, we can deduce that when rj > 1/2, the magnitude of the Casimir force is always a decreasing function of a. 

From ([5]), we can also deduce that in the high temperature limit, the Casimir force grows linearly with temperature, 
with leading term given by the sum of the terms with ^ = 0. In the low temperature limit, the Casimir force can be 
written as a sum of the zero temperature Casimir force and the temperature correction term, i.e., 



FSLsia) = F^^ia;T = 0) + ArF^^^ia), 
where the zero temperature Casimir force is given by 



(7) 



t-^ V ^4=2 I Li 

Fgi^{a;T = 0) = V cos{nk^r^) ^ '—K^ ( 2fcia 

TTa ^ — ' 



fci 



2a 



cos {nkii]) 



nki \ 



\ 



Kq 2kia^ 



d / , X 2 ' 
TTfci 



and the temperature correction is 



9 

TT a 



2aT 



(fc2,...,/cd)GN''-i 1 = 1 



fei=-oo (fc2,...,fcd)eN''-i / I 7r(fci--g) 




^■=2(1^)' 



exp 




- 1 



Observe that the temperature correction term goes to zero exponentially fast when the temperature approaches zero. 

The expression ([5|) for the Casimir force is suitable for studying the behavior of the Casimir force when a ^ Li, 
i = 2, . . . ,d and a ^ 1/T. On the other hand, the expression ((7|) is suitable for studying the behavior of the Casimir 
force when <C a <C 1/r, i = 2, . . . ,d. To study the behavior of the Casimir force when a <C 1/T <C Li or 
1/T <IC a <C Li, i = 2, . . . , d, we need alternative expressions. The derivation of these alternative expressions is a little 
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TABLE L The value of a;„ where B„{x) = 0. 



n 




n 




n 




1 


1/3=0.33333333 


2 


0.422649731 


3 


0.46165930 


4 


0.48067038 


5 


0.49023768 


6 


0.49508143 


7 


0.49752864 


8 


0.49876076 


9 


0.49937940 



bit tedious, but can be done using the same techniques as in [34] and 4(1|. We list these alternative expressions for 
the Casimir force in the Appendix. From (|A5|) . we find that when a ^ 1/T ^ Li, i = 2, . . . ,d, the leading terms of 
the Casimir force i^Cas(a) are given by 



cos(7rfci?7) , ^p,/^ 1 71" TT 



(8) 



The remaining terms goes to zero exponentially fast. The dominating term of the Casimir force when a <C 1/T <C Li 
is thus 



^d+l 



cos {■nkiif) 



fei=i '^i 



d+l 



(9) 



When 1/T <^ a ^ Li, i = 2, d, wc find from (|A8P that the leading terms of the Casimir force Fcas(a) are given 
by 



poo / N _ sr^ K-i^) 



E 



cos(7rfci77) 



aTE { ; — , . . . , — 

2 L2 Ld 



(10) 



The remaining terms goes to zero exponentially fast. The dominating term of the Casimir force when 1/T <§; a <C 
is thus 



aT 



2d-i 



fci=i 



cos(7r/ci77) 



(11) 



From ([9]) and pT|) . we find that the signs of the dominating terms of the Casimir force when a <C 1/T ^ T^ and 
when 1/T <C a <C ii are governed by the sign of the function 



COs{TTkx) 



[0,1], 



k=l 



when n = d+l and n — d respectively. When n is even, B„{x) is up to a constant, the Bernoulli polynomial of degree 
2n. For any integer n, the function Bn{x) has been studied in (34| . It was shown that this function is decreasing and 
has a unique zero in the interval [0, 1]. Moreover, this zero point is less that 1/2. The point Xn where Bn{x) = for 
n = 1, 2, . . . , 9 is given in TablelH One can observe that a;„ is an increasing function of n. Therefore, we can conclude 
that when Li ^ a, 1/T, i — 2, . . . ,d, and Xd < r] < Xd+i, increasing the temperature tends to change the Casimir 
force from attractive to repulsive. 

In FIG. [2] and FIG.[3l we show graphically the dependence of the Casimir force on various parameters such as a, T 
and r] when the dimension d is three. The figure on the left of FIG. [2] shows that when the fractional order 77 is equal 
to 0.47 and the temperature T is equal to 2K or 5K, the Casimir force can change from repulsive to attractive when 
we increase a. The figure on the right of FIG. [5] shows the dependence of the Casimir force on the fractional order 
77. In general, the Casimir force changes from attractive to repulsive when we increase 77. FIG. [3] shows that when 
rj = 0.47 and 77 = 0.48, increasing temperature can change the Casimir force from attractive to repulsive. 

As a conclusion, we have studied the finite temperature Casimir effect on a piston moving freely inside a rectangular 
cavity due to a massless fractional Klein-Gordon field, where Dirichlet boundary conditions are imposed on the walls 
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FIG. 2: Left: The Casimir force Fgi^{a) as a function of a when d = 3, L2 = is = 0.01m, t] = 0.47, and T = IK, 2K and 5K 
respectively. Right: The Casimir force F§^g{a) as a function of rj when d = 3, L2 = — 0.01m for various values of {a,T). 




T (K) T (K) 



FIG. 3: The Casimir force Fg-^^{a) as a function of T when d = 3, L2 = is = 0.01m, 77 = 0.47 (Left) and r? = 0.48 (Right), and 
a — 0.8mm, 0.9mm, 1mm and 1.2mm respectively. 



of the cavity and fractional Neumann boundary condition of order 77 is imposed on the piston. We show analytically 
that for T] larger than 1/2, the Casimir force is always repulsive. Moreover, the magnitude of the Casimir force is a 
decreasing function of a. This is another instance where we can obtain repulsive Casimir force in the piston setting 
besides the hybrid boundary conditions considered in j4l| and the Robin boundary conditions considered in 42] . The 
more interesting case is when rj is less than 1/2. We show that for the values of rj less than 1/2, the Casimir force 
can always become attractive when a is large enough. For certain values of rj less than but close to 1/2, it has been 
demonstrated graphically that the Casimir force can be attractive or repulsive depending on the aspect ratio of the 
cavity and the temperature. 



APPENDIX A 



First, we want to compute the A — > 0^ asymptotic expansion of the sum ^ up to the constant term in A. Let 
C(s; a) be the zeta function 



r(s) Jn 

fcGZxN'*-! /=-oo ^ ' ''^ fcGZxN'*-! /=-oo 



(Al) 
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Then 



d_ 

ds 



2s pc+ioo 



s=0 ^ J c—ioo 



^ r(z)A-^C 



z + 2s 



a dz. 



Using the formula 



ki— — oo 



A,'! — — oo 



(A2) 



one can show that the zeta function C{s',cl) (|Aip only has simple poles at s = i/2, 2 < i < d + 1, with residue 
Res,^.C(s;a) = "'t}T^^ , where 



E 

2<]l<j2<■■■<jd-l-^<d 



Therefore, C(0;a) = and for s e (-1/2,1/2), 

^/ r(z)A-c^^ 

From this, we find that 



d+l 



a\dz = 2Y, ^^^T^ri^^i^ - 2s)\^'-' + C(s; a) + 0(A), A ^ 0+. 



d+l 



^ ^ -2 ^ TSrrCd+i-^Ca) (log[AAi]2 - 2m) - C'(0; a). 

i=2 ^ \2) 



(A3) 



Next, we derive a formula that is needed for the computation of the Casimir force ([4]). Substitute (|A2p into (jAl[) . 
we find that 



C(s;a) = 



a r (s — i 

7^ r(s) ' 



/ 



4a 



^ ^ cos(7rfci?7) 



keN'' /=-ioo 



kia 



\ 



X -f'',-! 2fcia^ 



d / , X 2 
TTK," 



(27r/r)2 , 



where 



00 / d / , \ 2 



ci(^)= E E E(^) +(2-^^)^ 

(fe2,...,fcd)GN''-i ;=-oo \i=2 



Therefore, 



^C'(0;a) = -2aCi (-^) -4Rc ^ e'^'^'^" exp ( -2fcia. ; 



Trfc,: 



(A4) 



From this, we obtain the expression for the Casimir force ([5]). 

Finally, we list two alternative expressions for the Casimir force FQl^{a). To study the behavior of the Casimir 
force when a <^ 1/T <^ Li, i = 2, . . . ,d, we have 



(7r|fci-i|)^rV 

L 2 a 2 



tt/ fci — ^ I 



cos(7rfci??) a (-1)" ' ^2 2 ^ ^ 
2^ U+2 + 9d-l ^i+i '^J 



ki^l 1 j^O 



aT 



aTE 



1 TT 

2'I 



1 7- 1 ■ ■ • : 

2 -^d 



TT 
17 



d-1 



(-1)" 



-i-i 



J +2 



Cr{j + 2)T^"+2 
(A5) 



aT 
2^ 



d-i 



/ ^ .(fc,-a) 



/ 

2 

V 



2<ii<...<ij<c! 



fei = -oo (fcij ,...,ki.)eZ-i\{0} 



{2TTlTf 



27r/r 



xE E 



{2TilTf 



aT 
2^2 



E(-i)'-^~^' E 



2<ii<...<i,<d 



3 + 1 

TT 2 



(A6) 



where E [—^; . . . , j-j] is the analytic continuation of the function 



E s 



E 



,( 




2 


'irkd 






\ 


+ . . . + 







to s — —1/2, which is equal to the zeta regularized Casimir energy for massless scalar field inside a rectangular box 
[0, L2] X ... X [0, Ld] subject to Dirichlet boundary conditions (see e.g. [i^ 13, S^). It can be shown that (see e.g. 

Mm) 



El---— — 

2' L2 ' ' irf 



rf-1 v(i±^) 

i=i -^d 2<ii<...<ij-i 



■j-i 



^E(-i)'-" 

J=2 



E 



■i-i 



E E 



2TTkd 



4 L 



E (^»p-^»: 



(A7) 



To study the behavior of the Casimir force when 1/T <^ a ^ Li, i — 2, . . . , d, we have 



aT ^ {-ly-^-i jSj ^fj + l\ ^ cos(7rfci7y) aT ^ {-1^-^-^ . 



3=0 



TT 2 



00 00 



E 



1 = 1 ''l J=0 7^ 



E — HI — ■^'^j E E<=°^('''^i'?) 



fei=l Z=l 



/27r/ry 



if^ {iirhlTa) ~^Y1 ^^T+T— E E ^"^(^^i^) U {iiThlTa) 



d-i-j 



00 00 



{2TrlTy 



2d-A 



aT 
2^ 



X is:.- 



E("i)^-^-^ E 



E E 



(fcia)" 



2<ii<...<ij<£; 



(27rff)2 



V 



(27r?T)2 



aT£' 



1 TT 



2'L2"''Ld 2<^ 



2+1 



27r/r 



E _ _ _ 



\ p=i 



(A8) 
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